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Abstract 

We reexamine the long-standing problem of the microscopic derivation of a 
particle-core coupling model. We base our research on the Klein-Kerman ap- 
proach, as amended by Donau and Frauendorf. We describe the formalism to 
calculate energy spectra and transition strengths in some detail. We apply our 
formalism to the rotational nuclei ^^^'^^^Gd, where recent experimental data 
requires an explanation. We find no clear evidence of a need for Coriolis atten- 
uation. 



1 Introduction 

In this paper we use the Kerman-KIein method |[^, ^ Q to derive a microscopic core- 
particle model. This method is based on Heisenberg matrix mechanics, where exact 
eigenstates are used, but the matrix elements of operators are the unknown quanti- 
ties. In principle the formalism gives the same results as the Schrodinger equation. In 
practical approximations the method leads to interesting approximate models. Espe- 
cially when some matrix elements are known experimentally, we can construct models 
that are a hybrid between a phenomenological and microscopic model. In the case of 
core-particle coupling such an approach has been introduced by Donau and Frauen- 
dorf. We study the usefulness of this approach in a serious numerical apphcation to 
the rotational nuclei ^^^'^^''Gd, and extend the method to calculations of transition 
strengths. 

Our calculations, restricted to deformed nuclei in the current application, generalize 
the particle rotor model, and we shall compare our results to those of less microscopic 

*address after Sept. 1, 1993: Institiit fur theoretische Physik III, Universitat Erlangen-Niirnberg, 
D-91058 Erlangen, Germany 



1 



calculations using this model. The model, first introduced by Bohr and Mottelson |Q, 
has been used for the study of both low and high spin states of weakly and strongly 
deformed nuclei. For the past three decades it has continued to provide a framework for 
the analysis of the rotational band structure of odd nuclei. In spite of its approximate 
nature, it has been used successfully to explain both the structure of deformed and 
transitional nuclei 0. 

In its most primitive version one couples a single particle to a rotating (even-even) 
rigidly deformed core. The coupling of the single particle to the core is generally 
approximated by a deformed Nilsson potential. Due to the rotation there is also 
a Coriolis interaction between the total angular momentum of the system and the 
single-particle angular momentum. In later applications one finds extensions of the 
original model (we shall use calculations for ^^^'^^''Gd as references throughout the 
discussion). The Nilsson potential can be replaced by a more realistic deformed Woods- 
Saxon potential 0, pair correlations can be taken into account by a quasi-particle 
transformation of BCS type 0, ^ p, and mixing between different major shells 
(states with principal quantum number differing by two) |T0|, |ll|, |12|, |T3|, ||, |^ can 
be included. Finally one can allow for phonon excitations of the core |^, [1^, 15], in 
contrast with the standard core-particle coupling model (CPC model) where the core 
is assumed to be a rigid rotor. 

Even these extended models often suffer from an inadequacy, the so called Corio- 
lis attenuation problem. In order to reproduce the experimental energy spectra, the 
Coriolis matrix elements have to be reduced by up to 50%. This problem has been 
studied extensively over the years. Many physical effects previously omitted in the 
model have been invoked as the cure for this disease: the proper treatment of the two 
body recoil effect [0, the inclusion of octupole degrees of freedom ]T^, and the use 



of an angular- momentum dependent moment of inertia and pairing gap [18|. A more 



complete list of all possible suggested causes is given in Ref . [|I9[ . It is quite conceivable 
that one needs to combine several of these extensions to eliminate this problem. As 
yet, however, no clear consensus has emerged concerning the source of the attenuation 
problem. 

In the present paper we adopt the Kerman-Klein approach, emphasizing the steps 
required in the transition between a strict microscopic approach and one which contains 
compromises aimed at making applications easier to carry out. In effect we shall trace 
the steps of Donau and Frauendorf in their core-quasiparticle coupling model 



2T|. In their work a simple pairing and quadrupole-quadrupole separable interaction 



Hamiltonian was considered. Only a few application of this type exist, see Refs. p^ , 
and the review in Ref. |T9[. In our work we also introduce mixing between states with 
AA^ = ±2, where A^ is the principal quantum number. We believe that the pursuit 
of a microscopic approach will lead to a more complete CPC, and give some insight 
into the problems of phenomenological CPC models. In addition the microscopic basis 
allows for the straightforward inclusion of different effective interactions. 
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This paper is organized as follows. In Sec. 2, we discuss the basis of the model. 
First we describe briefly the basis of the Kerman-Klein method and then summarize 
the Donau-Frauendorf modifications |2y, ^|. We then turn to the application of 



a slightly extended form of the Donau-Frauendorf formulation to the nuclei ^^^'^^^Gd. 
These form an excellent testing ground for our model since extensive calculations have 
been performed for these nuclei. There also exists a body of recent experimental data 
on the transition strengths that has not yet been explored sufficiently Finally, 



in the last section, we give a review of the other theoretical descriptions available for 
these nuclei, and compare our results with these previous efforts. 

2 The Model 

2.1 Kerman-Klein Method 



The Kerman-Klein method was introduced in the early sixties [25| to provide a rota- 
tionally invariant description of deformed nuclei in the framework based on the equa- 
tions of motion (EOM). Among many applications, it was applied to a study of the 
foundations of core-particle coupling models Overall it has been applied to 



various nuclear many-body problems as well as to field theories [^, Here we 
use it as the starting point to derive and extend the CPC model. As has been stated 
in the introduction, one of the major difficulties with the CPC model is the Corio- 
lis attenuation problem. We believe that the Kerman-Klein method may be a good 
starting point for an investigation of this problem. The method starts from a shell 
model Hamiltonian and provides expressions for the properties of low-lying collective 
states of an odd-mass nucleus and its even-even neighbors. These expressions include 
the energy eigenvalues, the single-particle coefficients of fractional parentage (cfp) and 
matrix elements of the electromagnetic operators. In practice, Hamiltonians with only 
simple ingredients have been studied (in the simplest case single particle Hamiltonian 
plus monopole pairing and quadrupole-quadrupole interactions), and only relatively 
few low-lying collective states in the even cores have been included (in the simplest 
case only the ground-state rotational band). Since the Hamiltonian must include ex- 
pressions for both the long-range part and the short-range part of the nuclear force, 
we start with a Hamiltonian containing multipole expansions involving particle-hole 
as well as particle-particle contributions. This separation simplifies the application of 
the EOM method. 
We thus write 



H = Y^ haaltto, + ^ H H Facdb{L) B[j^j^ (a, c)Blml {d, h) 

a ° abed LMl 
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+ ^ H H Gabcd{L) Al^^{a, b)ALMAc, d). 



abed LMl 



Here ha (a = (ja, ^a) and a = ja) are the spherical single particle energies, the operator 
Blm is the particle- hole multipole operator, 



al X Qb 



mam,, 



Ml ' 



and Alm is the particle-particle multipole operator 



mami. 



L 

Ml 



(2) 



(3) 



where are the Clebsch-Gordan coefficients and 
F are the particle-hole matrix elements 

Facdb{L) = 2^S^Sp C^^'^ ^5^^ Vafi^yS 
m's 

and G the particle-particle matrix elements 

^abcd[^) = ^ Z^^a/S ^a/3tS- 



{—Y" The coefficients 

(4) 

(5) 



The task we set ourselves is to obtain equations for the states and energies of the 
odd nucleus assuming that the properties of the even nuclei are known. The states 
of the odd nucleus (particle number N) are designated as li^J/j) where u denotes 
all quantum numbers beside the angular momentum J and its projection /i. The 
eigenstates and eigenvalues of the neighboring even nuclei with particle numbers (A^±l) 
are \IM n {N ± 1)) and E^^^, respectively, where n plays the same role for even nuclei 
as u does for the odd nuclei. The equations of motion (EOM) are obtained by forming 
commutators between the Hamiltonian and single fermion operators. 



[aa,H]^haaa + jY^Y^'^a^^ C;acbd{L)a\ALM{c, d) 



+ 7 X! X] SjC^^ Facdb{L)a-yBLM{d, b), 

^ bd-y LM 



(6) 



^1 = -haai 



YYCai'Gacbd{L)A[^{c,d)a, 



bdj LM 



7 E E s,C^^Facdb{L)Bl^{d, b)a\. 



(7) 



bdj LM 
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The matrix elements of these equations provide expressions for the single particle 
coefficients of fractional parentage U and V, defined as 



?7}^(a; IMn) = (z/J/x | at \ IMn{N - \] 

IMn) = {uJfi \ a„ | IMn{N + 1)) . (9) 



These are obtained by the evaluation of matrix elements of the EOM (^0) between 
eigenstates \yJ fj) and \IMn{N ± 1)) and the use of the completeness relation (details 



of the derivation can be found in |^). The resulting equations can be cast into the 
block matrix form, 




-t + uo 





(10) 



Here e are the single particle energies measured relative to the chemical potential, A, 
the elements of uj are the excitation energies in the even nuclei, and F and A are pair 
and multipole fields. More specifically we have 



€aIMn,"/I'M'n' — {^N — ha)5a-y5 1 Ii5 M M' ^nn' ■, (H 
1 



= ^[Er'-K-'\^ (12) 

^aIMn,-/I'M'n' = i^In ^ ~ ^0 ^)^o,-yhl'5 MM' ^nn' i (13) 



^ aINIn,-iI'M'n' ~ X! ^■y^a^ ^ FgcdbjL) 
LMl bd 

X {I'M'n'{N ± 1) I BlmM b) I IMn{N ± 1)) , (14) 

^aIMn,jrM'n — X! ^ G acdbiL) 
LMl bd 

X {I'M'n'{N - 1) I AlmA^, b) \ IMn{N + 1)) . (15) 



In order to specify the solutions uniquely a normalization condition is required. 
From the anti-commutation relations of the fermion operators, one takes the appropri- 
ate matrix elements using a sum over intermediate states and finds 



^ E [|f/;,(«/M;n)p + |V-(a/M;n)p] =1, (16) 

^ ' a IMn 

where 

r] = ^(2j, + i). (17) 

jo. 



5 



All of the above equations are still exact and represent a formulation of Heisenberg's 
matrix mechanics. The fundamental practical problem is to find a suitable 'collective' 
subspace of the complete Hilbert space. We assume that this 'collective' subspace 
decouples from the remaining states, at least approximately. This requires that any 
matrix element of the relevant multipole operators between a collective and a non- 
collective state is negligible. 



2.2 The Donau-Frauendorf Model 

In this section we introduce the approximations made by Donau and Frauendorf, and 
use these as the practical basis for our investigation. First we assume that we have a 
finite sum of separable interactions, where each of the terms in and is of the 
form 

FacdbiL) = -2KLFac{L)Fdb{L), (18) 

GacdbiL) = -2gLGac{L)GM. (19) 

Here k^, gi are the strengths of the corresponding multipole and pairing forces, and 
Fac{L) and Gac{L) are the appropriate single particle matrix elements. For the moment 
we consider only monopole-pairing and quadrupole-quadrupole interaction (the inclu- 
sion of hexadecapole and the quadrupole-pairing interactions is left for later work). 
The quadrupole-quadrupole and monopole-pairing operators are 

QM = E^-(2)52M(a,c), (20) 

ac 

A = E^-(0)Aoo(a,c). (21) 

ac 

co-workers have adopted Next we identify the "collective" subspace appropriate for 
the applications that follow. The obvious choice is the ground state rotational band 
\IMK = 0), where K is the projection of angular momentum on the body fixed axis. 
The fact that this "collective" subspace approximately decouples from the remaining 
states is equivalent to the requirement that inter-band BE{2) transitions are much 
smaller than the intra-band transitions. This is generally true for deformed nuclei. A 
last approximation for deformed nuclei, convenient but not necessary, is to replace the 
excitation energies as well as the multipole and pairing fields by their average values 
for the two neighboring even-even nuclei. 

Before we go on to the next step it is convenient to rewrite the EOM with reduced 
matrix elements, using the Wigner-Eckart theorem. We can thus derive equations of 
motion, again dropping all indices. 



/ e + u + T A 

A -e + uj-T 




( ^ 

u 



(22) 
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where V and U are the reduced cfp 

f/)(l;/)a = (z. J II at II /(AT - 1) 

Vj\l;I)a = (z/J||Sj|/(iV+l)), (23) 
the self-consistent fields are represented by the expressions 

^ai,ci' = -l^2{-y^^'^'l^ ^; J } (/ ii Q2 11 /') fu2), (24) 

and the generalized pairing force is approximated by a constant gap energy for all 
levels, 

Aai,ci' = (/|| A||/') 5ac ^ (Oil A||0) 5n'5ac. (25) 
The normalization condition becomes 



W[\UXl;I')a\'+\Vj\l;l')a\']=l. (26) 



al 

Here the inclusion of only the collective subspace is somewhat more questionable than 
for the EOM, but there appears to be no simple alternative. 

The final step in the derivation is the replacement of the self consistent fields by 
phenomenological inputs. In this context (/||Q||/') and A are the quadrupole matrix 
elements of the neighboring even nuclei and A their gap energies. The quadrupole 
fields can either be extracted from experiment or calculated using a phenomenological 
model such as the Bohr-Mottelson model IHI, 



(/'iiQii/) = go^Y6^V(2^ + i)(2/' + i)(^;;;j, (27) 

where go is a phenomenological input which can be obtained from experimental values. 

Having described the assumptions involved in the derivation of the model, we turn 
to the task of solving the resulting equations. The main difficulty in solving those 
equations is that the set of solutions ( ^2]) is over-complete by a factor of two. This 
is a consequence of the fact that the basis states \I{N — 1)) and a \I{N + 1)) form 
an over-complete and non-orthogonal set as in the standard BCS model. Half of the 



states found by solving (p^ are spurious, and these have to be identified. To remove 



the spurious solutions one can proceed in the following way |^, First build the 
density matrix 

uJ 

{I{N - l)\K\\uJ) (z/J||al,||/'(iV - 1)) (/(iV + l)H\\uj) (z/J||al,||/'(iV - 1)) 
(/(iV-l)||a„||z/J) (z/J||a,,||/'(iV+l)) {l{N + l)\\al\\uj) {uJ\\aa'\\r{N + 1)) 



:28) 
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which can also be expressed in terms of the reduced cfp U and V, 

U^l; I)aU^il; I')a' f/}(l; I)aV^{l; I')a' 

/)af/}(l; /')«' V}{1- I)aV}{l; I')a' 
This density matrix has two important properties; it commutes with Hamihonian 



J Jv 



■ (29) 



[H,D] 



0, 



(30) 



where H is the matrix Hamihonian of Eq. (^) and the scaled quantity il D acts like 
a projection operator 

{nnf = nD. (3i) 

(Here (^) follows from the EOM and (|3l|) expresses the ortho-normalization conditions 
that follow from (^) and (gB|).) 

Therefore, D has eigenvalues 1/Q or 0. The eigenstates of the Hamiltonian which 
are also eigenstates of D with eigenvalue 1/Q describe the real quasiparticle excita- 
tions in the odd-A^ nucleus, whereas the eigenstates with eigenvalue characterize the 
spurious states. 

We next decompose the Hamiltonian as. 



H = Hqp + H, 



where 



qp 



-e-r A 
A e + r 



UJc 

oj. 



(32) 



(33) 



Hqp is interpreted as a generalized quasi-particle Hamiltonian and as the core Hamil- 
tonian. H is separated into two parts because Hqp is antisymmetric with respect to 
particle- hole conjugation (aj — Oq ), and therefore, as in the usual BCS theory [^], the 
solutions of this part of the Hamiltonian are divided into two sets with sign-reversed en- 
ergies; the negative energy solutions are the spurious states. This decomposition gives 
a good starting point for a numerical technique to remove the spurious states. First, 
we turn off the core Hamiltonian He- We solve the eigenvalue problem, Hqp^^ = £°\E'°, 
and find the solutions \E'j. and where 



From the physical states \E'^ we can construct the density matrix D'^ as in 
the core Hamiltonian is gradually turned on. 



. Then 



H{j)=Hqp + jH, 



(34) 
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where 7 is a scaling parameter. The next step is to diagonahze H{'y) and find the 
new solutions, ^(7). If the change of 7 is not too big, the quantity '^^{j)D'^'if{'y) will 
either be close to or to l/fi. When the value of ^''"(7)D°^(7) is close to l/fi it is 
a real state and when it is close to zero it is a spurious state. Now we can construct 
the new density matrix from the physical states ^'+(7) which we call D^. Then 7 
is incrementally increased and this process is repeated until 7 = 1. In practice this 
method works very well since the change in the eigenvalues of D never exceeds 5%, 
even with only 15 incremental steps. 

When 7^0, which is the limit that the moment of inertia X goes to infinity, 
the core excitations vanish (adiabatic limit). In this limit, Hqp is equivalent to the 
usual Nilsson plus pairing Hamiltonian. To demonstrate this, we first diagonalize the 
particle- hole part (e + F) of H^p. The resulting eigenvalues — A are characterized 
by the index k. For each different k, F defines a different fixed quadrupole field or, in 
other words, an instantaneous orientation of the core. The index k can be related to 
K, the projection of the angular momentum on the body fixed frame, by a method 
described below. Finally, the solutions of the Hqp, including pairing, are 



The remaining problem is to find the correspondence between the values of K and 
those of k. This is not straightforward because Rgp is an angular momentum conserving 
quantity; therefore K can not be obviously associated with the eigenvalues. We use the 
fact that among states with different J, states with the same value of K have identical 
eigenvalues. We can extract a unique state of maximal K value in the following way: 
Compare the J + 1/2 doubly degenerate solutions for angular momentum J with the 
J — 1/2 doubly degenerate solutions for J — 1. Since the maximum values of K are J 
and J —1, respectively, the additional pair of solutions for angular momentum J must 
have \K\ — J. We then proceed with J — 2, and so on. 

2.3 Electromagnetic Transitions 

The electromagnetic multipole operators are one-body matrix elements of the form 



where e^-y is either the effective charge or the effective g factor, and ta-y are the single- 
particle matrix elements of the transition operator. 

It is convenient to decompose the electromagnetic observables into core and particle 
contribution, as in the conventional CPC. To do this we must express the matrix 
elements of the transition multipole operator of the odd nucleus, {J'ij,'i''\Tim\Jij>i'), in 
terms of the matrix elements of the even nuclei (which are known from experiment) 




(35) 




(36) 
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and the known eigenfunctions U and V. We use the decomposition of the eigenstates 
of the odd nucleus |20| , 



iJf'^) = ^ E [UjM, IM) at \IM) + V^ia- IM) 



IMa 



IM 



(37) 



where 



\IM) = \IM{N-1)), 
Tm) = \IM{N+1)). 



(38) 

We then rearrange the order of the single particle operators a and and the collective 
operator T using their commutation relations and utilize sums over intermediate states 
of the even nuclei such that operator T stands between states of the even n uclei and 
the single particle operators occur between an even nucleus state and an odd nucleus 
state. Finally we have 



{J'fi'i^'inmlJH = ^ E E [f/j^(«;/M)f/5,;,(a; I'M') (I'M' I TLm \ IM 

IMa I'M' 

+V^a; IM)Vj:^,{a; I'M') li'M' \ T^m \ IM 



1 



(39) 



+ ^ E E ^.7^-7 iMW;^,{r, IM) + viia- iM)Vj%,ir, im) 



IMa 7 



Equation ( pQ] ) clearly exhibits a separation into particle and core contributions. With 
the definition of the reduced matrix elements (^) we ultimately obtain 



[j'u'WTaju) 



+ 



n 



E innwi') 



air 



E( 

lac 



\ja + J' + I 



ya+J+i' 

' J J' 2 

ic ja I 



r I 2 

J J' ja 



U:;{a; I)U'j'{a- 1') + V7(a; I)V^',{a- 1') 



Uj{a- I)U:;:{c; I) + Vj{a- I)Vj: {c; I) 



where we assume that the matrix elements of Tl do not depend on A^, 



{I\\Tl\\I') 



I\\Tl\\I' 



(41) 



For the electric quadrupole transitions studied in this paper, the matrix elements 
(/||Tl||/') become the matrix elements of the quadrupole operator The single 

particle matrix element tac{L) are 



Fac{L) 



a II rV^ II c 



(42) 



For these electric transitions we have only one parameter to fix, the effective charge. 
As we will see later this is a minor problem since the single particle contribution is 



(40) 
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much smaller than the core contribution. The quadrupole matrix elements of the core 
can be either derived from the Bohr-Mottelson model as in ( ^7]) or extracted from the 
experimental results. 

For the magnetic (Ml) transitions the matrix elements of the operator is 



{I\\Tl\\I') = {I\\M\\I) deltajp, 
and the single particle matrix elements tac become 

tac = Gac = (ja II M II jc) 



An 



ly 



x{l-k) 



1 , k., 

-gs-gi{l + -^ 



with 



k = Ua + l)(-l)(>+'»+l/2) + a + i)(-l)(>+'=+l/^). 



(43) 



3(2j, + l)(2j, + l), .^,W/2/^ Ja 1 Jc 



-1/2 1/2 



(44) 



(45) 



3 Calculations for 



155,157 



Gd 



We decided to concentrate the calculations on the nuclei ^^^'^^'''Gd, since recent detailed 
studies of the transitions seem to require explanation. The even cores i54,i56,i58q^ 



are well deformed and have almost rigid rotational bands. Therefore, the quadrupole 
field r can be estimated using the Bohr-Mottelson model as 



22.21 



al,cl' 



' ] I I' J 



^(2/+l)(2/' + l) 



1 r 



2 



(46) 



where qac are the single particle quadrupole matrix elements (labeled before by Fac{2)) 
in fm. The factor 22.21 comes from the transformation from the intrinsic quadrupole 
moment to the deformation parameter p. The core energies are loj = where for 

the moments of inertia the values X^^^ = 0.0122 MeV , J^^"^ = 0.0124 MeV are chosen. 



For the pairing potential we used the value A = MeV. The chemical potential A 



was calculated according to (|1^) from the difference between the ground state energies 
of the cores. The values found are A^^^ = -7.4845 MeV and A^^^ = -7.1488 MeV. The 



experimental values for the quadrupole deformations of the cores are |30 



154 


Gd : 


/520 


= 0.3294, 


156 


Gd : 


/320 


= 0.3378, 


158 


Gd : 


020 


= 0.3484. 
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Figure 1: Single particle eigenvalues for Woods-Saxon potential . 



For the odd nuclei we used the averages of the two neighboring even nuclei, (3\q^ = 
0.329, Plo'^ = 0.3431. The strength of the quadrupole force k was taken as a free pa- 
rameter. The single particle wave functions were found from a Woods-Saxon potential 



-V. 



1 + exp 



(47) 



including spin-orbit interaction with parameters a = 0.65, Rq = 1.25^4^/^ and Vq = 
51/ (l — l ^p^^ ^ MeV. The single particle quadrupole matrix elements qac take values 



according to the standard formula 



Oa II Q II Jc) 



/5 (2j,-t-l) (2j, 

4:71 



^(_l)ic-l/2j 



ja 2 jc 

-1/2 1/2 



(48) 
(49) 



The resulting single-particle levels are shown in Fig. |I[ 

First we apply the method to positive parity states of ^^''Gd. The experimental 
results identify the 5/2^ [642] band as the lowest positive parity band, followed by the 
3/2+ [402] and 1/2+ [400] bands at excitation energies ~ 0.4 MeV and ~ 0.6 MeV r 
espectively. Figure (||a) shows the solution in the adiabatic limit (7 = in Eq. (0)) 
as a function of deformation using only the 1^13/2 quasi-particle state. As stated before. 
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Single Particle States Rotational states 




Figure 2: Energy levels for the positive parity states of ^^''Gd. Only single-particle 
states from = 6 major shell are included. The 5/2 [642] is below the 3/2 [402] and 
l/2[400] at deformations larger than 0.33. The quadrupole-quadrupole strength is set 
at K = 0.311 MeV/fm. The black squares represent the experimental values. 

this should be equivalent to the Nilsson diagram (including pairing). At zero deforma- 
tion we have a degeneracy because of the spherical symmetry. As the deformation is 
increased the spherical symmetry is broken and the degeneracy is lifted. As a result 
there are 14 states each one with different K, taking values from j to —j. Because of 
the signature symmetry, states with K having opposite signs are degenerate. For the 
case of prolate deformation, quasi-particle states with smaller K are lower in energy. 
In the case of quasi-hole states the opposite is true. We found that for quadrupole- 
quadrupole strength k = 0.31 MeV/fm^ the band-head energies are in best agreement 
with experiment. The 5/2+[642] state is below the 3/2+[402] by about 0.4 MeV, but 
the 1/2+ [400] is off by 1 MeV. In Fig. we show the result when the core excitations 
were turned on. The lowest band 5/2"*" [642] band reproduces the experimental band 
and to the very limited extend that data is available, so does the 3/2"*" [402] band. The 
1/2+ [400] band is off by 1 MeV. 

For better results we include more single particle and single hole states ; two more 
major shells were included, the = 5 and N = 7 (see Fig. |l]). For k = 0.23 MeV/fm^ 
the band-head energies are in good agreement with experimental values. This change 
in the value of k was to be expected, since we are using an effective interaction. It is 
therefore natural to need different strengths for the interaction at different dimensions 
of the included single-particle space. The zig-zag shape of the 3/2+[402] and l/2+[400] 
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Figure 3: The positive parity energy levels of ^^^Gd. The circles correspond to the 
theoretical predictions and the squares to the experimental values. The strength of 
quadrupole-quadrupole interaction is k = 0.23 MeV/fm^. Here we include states from 
3 major shells. 

rotational bands reveals that the Coriolis interaction is strong compared to single parti- 
cle excitations (quadrupole-quadrupole interaction) such that it creates the staggering 
effect. Since there are no experimental data for those bands, we cannot draw any firm 
conclusions concerning this behavior. The fact that the 1/2+ [400] band- head is still off 
by 0.3 MeV may imply that a more sophisticated interaction is needed to fully describe 
the structure of this band. 

The next step is the application of the method to the negative parity states where 
more data is available. The experimental situation can be summarized as follows: 
The 3/2~ [521] is the ground state band. Then follows the ll/2~ [505] hole band and 
the 5/2" [523] particle band. The 3/2" [532] band and the 1/2" [530] band, which 
are not certain, are higher in energy and almost degenerate. If we only include the 

= 6 major shell, the strength k should be set at a value of about 0.2 MeV/fm^ 
for the 3/2~[521] band-head to be the ground state, which agrees with all previous 
calculations [^ Figure ||a shows our "Nilsson" diagram. Again particle states 
with lowest K have lowest energy. The states deriving from the hole states obey 
the opposite rule; the lowest state has the highest K value. In addition to this, we 
note that quasi-particle and quasi-hole states originating from different single particle 
orbitals interact weakly with each other. This leads to narrow avoided crossing between 
states with the same K value. States with different principle quantum number N also 
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Figure 4: Negative parity states of ^^^Gd. The squares correspond to experimental 
values. Only states from the N = Q major shell are included. The strength of the 
quadrupole force k is set to 0.201 MeV/fm^. 



interact weakly. This is because the matrix elements (r^) for states having AA^ not 
equal to zero is smaller by a factor of the order of 10 as compared to states having 
the same N. Quasi-particle (or quasi-hole) states having angular momentum differing 
by four units and having the same K couple strongly. For Aj equal to 2 or their 
interaction is weak. This is a consequence of the relevant geometric factors (three and 
six-j symbols). For example, the quasi-particle states of /5/2 repel the hQ/2 states with 
the same K; as a result the band-head energies are quite sensitive to the position of 
the single-particle energy of the /5/2 orbital. Because of the approximate nature of the 
Woods-Saxon potential, some adjustments to the single particle energy levels will be 
necessary for better results. The 11/2^ [505] single hole state deriving from I/111/2 does 
not interact with any other state. We can easily force the 11/2^ [505] band- head to be 
0.4 MeV above the ground state by lowering the l/iii/2 single-hole energy level by 1.5 
MeV. Figure shows the calculated band structure compared to the experimental 
values. The ground state band 3/2"[521] and ll/2~[505] band-head are reproduced 
very accurately. The 1/2" [530] band-head is in good agreement with the experimental 
value but the rotational structure is distorted. This is due to the strong Coriolis 
interaction. The 5/2" [523] and 3/2" [532] band-heads ares way off (~5 MeV and ~10 
MeV respectively). At this point we can no longer improve the prediction by small 
adjustments of the single-particle or single-hole energies. 

In an attempt to improve agreement with experiment we include all single particle 
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and hole 1 evels from from the = 4 major shell as well as intruder states from 
N=5 and N=7 shells. As was to be expected, the strength of the interaction must be 
adjusted because of the change in the size of the single particle space. For best results 
the quadrupole-quadrupole strength was found to be k = 0.397 MeV/fm^. Though it 
is clear why /t is different for different dimensions, it is not clear why k is not the same 
for positive and negative parity states. For one major shell the positive parity states 
require k = 0.31 MeV/fm^ whereas for the negative parity states k, = 0.20 MeV/fm^. 
For three major shells the positive parity states require k = 0.2 MeV/fm^3, whereas 
for the negative parity states k = 0.39 MeV/fm^. For the moment, this remains a 
puzzle. 

The results for this more complete calculation are shown in Fig. § and are more 
satisfactory than those for one shell. The 5/2~[523] band-head is at the right energy 
and the same is true for the l/2~[530] state. The ll/2~[505] state is at the right 
position without any change to the single-particle energy spherical energy found by 
solving the Woods-Saxon potential. Only the 3/2~[532] is off by ~ 1 MeV. The only 
problem is that the structure of the rotational band l/2~[530] is slightly distorted. 
This again may indicate that higher order interactions are needed to describe states 
of higher energy. The results are very satisfactory keeping in mind that very few 
free parameters were used. The most important result is that we can reproduce the 
rotational structure with accuracy better than any previous work, without using any 
attenuation or any other forms of interaction. At the same time the inclusion of the 
latter remains a relatively straightforward possibility 

In Fig. ^we show a similar calculation for the negative parity levels of ^^^Gd. Single- 
particle states from the A^ = 4, A^ = 6 as well as intruder states from the N = 5 and 
N = 7 shells were included. The experimental structure of the bands is very similar 
to those of ^^''Gd with the slight change in the band-head energies. The quadrupole- 
quadrupole strength was found to be k = 0.377 MeV/fm^, which is comparable to the 
one used for the negative parity states of ^^^Gd. The 5/2~[523] and the ll/2~[505] 
band- heads are reproduced well. The same is true for l/2~[530] but the rotational 
structure of the band is distorted. As for ^^''Gd, the 3/2~[532] is off by ~ 1.5 MeV. 
The results are as good as for ^^^Gd and the same conclusions as above apply for this 
nucleus. 

We have also calculated the electromagnetic transitions and compared them to ex- 
periment. In Fig. we show the experimental BE{2)'s for ^^^Gd and ^^^Gd compared 
to the theoretical calculations using (|27|). The theoretical calculation represents the 
average of ^^^Gd and ^^^Gd. As can be seen from the figure the model reproduces the 
experiment reasonably well with some evidence of deviations at the largest angular 
momentum. For odd nucleus states having angular momentum J < 19/2, the contri- 
bution of states in the even cores with angular momentum / > 12 will be small. In 
the case of ^^^Gd and ^^^Gd we have a similar situation. Therefore, we can use either 
the theoretical calculations or the experimental values. Various aspects of the BE {2) 
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Figure 5: Negative parity energy levels ^^^Gd. Single particle and hole states from 
three major shells are included. Squares correspond to the experimental values and 
circles to theoretical. The strength k, of the quadrupole force is 0.397 MeV/fm^. 
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Figure 6: Negative parity energy levels for ^^Gd. Single particle and hole states 
from three major shells are included. Squares correspond to the experimental values 
and circles to the theoretical calculations. The strength k of the quadrupole force is 
0.377 MeV/fml 
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Figure 7: The experimental values of BE{2)^s in two of the even-even Gd isotopes 
compared to the values calculated using Bohr-Mottelson model. 

for transitions from J to J — 2 and from J to J — 1 are illustrated in Figs, 



The experimental data are taken from p4|. The salient fact is (as was expected) that 



the core contribution is much larger than the single particle contribution (compare 
the scales in Fig. ^ with those in Fig. We see that the numbers of single-particle 
states included only affect the single-particle contribution but not the core contribu- 
tion. This shows that we are close to the extreme strong coupling limit, where the 
individual valence neutrons do not alter the collective behavior (core behavior). On 
the other hand the single particle contribution is sensitive to the number of single 
particle states, but if more than the major shell is included, this contribution to the 
BE{2)^s converges (Fig. We have also calculated the BE{2)^s for different values of 
the effective charge. The best fit is for close to one. However, since the contribution 
of the single particle term is very small, the value of Cefr is of little importance. The 
final results of theoretical calculations for ^^^Gd and ^^^Gd are presented in figure |^ 
and |l^ together with experimental data and results of other theoretical models. 

We next calculate the magnetic transitions. Here we have four parameters to fit. 
Because of the limitation on available experimental data for the magnetic moments of 
these states of ^^^'^^^Gd, we have to use a phenomenological model for their values, |^ 

{I\\M\\I)=gJ + g2l', (50) 

where gi and g2 are free parameters which we fit to the experimental points. In practice 
the experimental values permit a certain flexibility in the choice of these parameters. 
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quantity theoretical fit ^^^Gd ^^''Gd 

gl 0.28 ±0.05 027 OX" 

c/2 0.018 ±0.008 0.019 0.02 

gL 

qs -3.826 -3.586 -3.576 



Table 1: The values of the parameters used for the calculations of M(l)'s (last two 
columns). The second column shows the theoretical values for ql and gs, whereas the 
third column shows the results of a fit to ^^^Gd for the parameters gi and g2- 

The values found by fitting the experiment in ^^^Gd and the values used are shown 
in Table ([1|). We further adjusted gi, g2 and gs slightly in order to find the best 
reproduction of the experiment (see Tab. (|l|)). Figure [Tl| shows the calculations for 
the core contribution, the single particle contribution and the total Ml. As can be seen 
from the figure, the core and single particle Mi's have roughly the same amplitude, 
and their amplitudes where chosen so as to achieve the best total Ml. 



4 Review of previous Work 

We turn now to previous theoretical work on ^^^'^^''Gd, in order to summarize it and 
compare it with our own work and with experiment. We can divide this work into 
three major groups according to the formulation applied. 

In the papers of the first group, Refs. ||^, |^, |15[ only band-heads are calculated, 
using the quasiparticle-phonon coupling scheme. In this work, the underlying physics 
is that of the standard core-particle coupling (CPC) model. Allowance is made for 
phonon excitations of the core, which can have both a quadrupole and hexadecapole 
deformation. The single particle states are found using a Woods-Saxon potential, and 
one uses a BCS approximation to include pairing correlations. In addition, interactions 
between states with principle quantum number differing by two are also included. 
Only band-heads have been calculated with this model, and the results show little 
improvement, as regards experiment, over calculations based on the Nilsson model, 
although fewer parameters are used in this approach [|l^]. In all these calculations the 



Coriolis interaction is neglected. The argument is that the Coriolis interaction has 
only a slight effect on the energies and structures of the lowest parts of the rotational 
bands |7[. The results of the best two calculations of this type are compared with 



experiment and with our calculations in Fig. 12. 



The second group, Refs. [[L^ , p^ , |, || contains complete calculations of the rota- 
tional spectrum in the standard core-particle coupling model. Early works based on 
this model utilize a rotating deformed core which is coupled to a single odd nucleon. 
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Figure 8: The single particle contribution to the BE{2) ^^''Gd for different numbers 
of single particle states, using CeS = 1- 
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Figure 9: Comparison of our results with other models for ^^''Gd. The formulas for 
the Particle- Rotor Model were taken from Ref. M, for the Cranking model from Ref. 



2J]. The points with error bars are the experimental data p4[, the short-dashed line 



is the Particle-Rotor model, the long-dashed line is the Cranking model and the solid 
line is the present work. 
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Figure 10: Comparison of our results with other models ^^^Gd. See Fig. 
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Figure 11: The calculated values for the Ml for ^^^'^^^Gd. The parameters used are 
are listed in Table 1. The points with error bars are the experimental data P^|, the 
short-dashed line is the single particle contribution, the long-dashed line is the core 
contribution and the solid line is the total calculated M(l). 
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Figure 12: Negative parity band-head energies for ^^"^Gd from: (a) experiment, (b) 
reference 0, (c) reference (d) our work. 



The mean field from the core is represented by a Nilsson potential. The coupling 
between the core and the odd nucleon, the Coriolis coupling, is treated as a pertur- 
bation (strong coupling limit). In Ref |13] a deformed Woods-Saxon potential is 



considered. Again, pairing is taken into account by means of a BCS transformation. 
Another extension of the original model takes into account mixing between different 
major shells [|I2], |13[. In Fig. ^ we show the experimental values compared to three 

(Fig. ^3|d) is from [12|. In this paper the model 



theoretical models. The first one (Fig. pTSp) is from 
Hamiltonian includes a quadrupole and hexadecapole deformed core. The pairing is 
treated in BCS approximation. The deformation parameters f32o, Pio and the moment 
of inertia X are taken as freely adjustable parameters. The set of results shown in 
Fig. |l^c is from Ref. |T^. In this paper the same Woods-Saxon potential and BCS 



approximation are used. The model contains two additional free parameters: the mo- 
ment of inertia and the Coriolis attenuation factor k, which was found to be 0.63. The 



last part of Fig. ^ is based on our results, without any Coriolis attenuation. 

The third group includes models based on dynamical symmetry arguments; the 
pseudo SU(3) ||3^ and the interacting boson-fermion approximation (IBFA) 0]. The 
results are reasonable for both models but they will not be reviewed because they 
describe only the lowest bands. 
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Figure 13: Negative parity states for ^^''Gd from: 
(c) reference fl^, (d) our work. 



(a) experiments, (b) reference |T2[, 



5 Conclusions 

In this paper we have apphed to ^^^'^^^Gd the Core-Particle model as originally for- 
mulated by Klein and co-workers and modified for practical application by Donau and 
Frauendorf. When levels from several major shells were included, the band-heads were 
predicted quite satisfactorily. The rotational structure is predicted with high accuracy 
and only the higher bands are distorted. We believe that the reason for this distortion 
is the Coriolis interaction which is included naturally in the core Hamiltonian {He). 
The standard approach to this difficulty is to attenuate the Coriolis interaction |jT2|, |13| . 
Since there is no justification for this ad- hoc procedure, the best approach is to include 
other forms of interaction which will counteract this effect. The suggested candidates 
are the quadrupole-pairing and hexadecapole-hexadecapole interactions, which will be 
investigated in the further development of this work. 

The electric transitions BE (2) were predicted and the results are closer to the 
experimental values than any other theory for which comparisons have been made. For 
the magnetic transitions M(l) we reproduce the experiment but need 4 parameters to 
do so. These parameters are not fully constrained by experiment. 

Finally, it may be useful to summarize in general what has been accomplished and 
how the model may be extended. In our approach, the kinematical Coriolis coupling 
between the particle and the core is automatically included. There is in fact no natu- 
ral way of introducing an ad-hoc attenuation, for which we have found little evidence. 
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There is some indication in higher excited bands, but not enough data is available 
to substantiate this. In our calculations an excessive influence of this coupling can 
be dealt with only by modifying the effective interaction. In addition non-adiabatic 
corrections to the treatment of the core may be studied, either by use of experimen- 
tal matrix elements or by calculating higher order corrections to the Bohr-Mottelson 
theory. Finally to study experimental results at higher angular momentum, we must 
include excited bands of the core. Because of the modest success of the calculations 
reported in this paper, we feel encouraged to continue our investigations along the lines 
suggested above. 
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